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Phenomena associated with the controlled collapse of 
a Bose-Einstein condensate described in the experiment 
of Donley et al '1 / are explained here as a consequence 
of the squeezing and amplification of quantum fluctua- 
tions above the condensate by the condensate dynam- 
ics. In analyzing the changing amplitude and particle 
contents of these excitations, our simple physical pic- 
ture provides excellent quantitative fits with experimen- 
tal data on the scaling behavior of the collapse time and 
the amount of particles emitted in the jets. 

In the experiment described by Donley et al. pQ, a 
Bose-Einstein condensate (BEC) in a cold (3nK) gas 
of Rubidium atoms is rendered unstable by a sudden 
inversion of the sign of the interaction between atoms. 
This is done by altering the binding energy at Fesh- 
bach resonance with an external magnetic field. Af- 
ter a waiting time t co u apse , the condensate implodes 
(called Bose-Nova), and a fraction of the condensate 
atoms are seen to oscillate within the magnetic trap 
which contains the gas. These atoms are said to be- 
long to a 'burst'. After a time r evo i ve the interaction 
is suddenly turned off. For a certain range of values 
of revolve i new emissions of atoms from the condensate 
are observed. They are called 'jets'. Jets are distinct 
from bursts: they are colder, weaker, and have a char- 
acteristic disk-like shape. 1 



* Emails: calzetta@df.uba.ar, hub@physics.umd.edu 
1 We call attention to the distinction between the 'Bose-Nova' 
[Tl experiment studied here and other BEC collapse experiments 
At magnetic fields around 160G, where the effective scat- 
tering length is of the order of 500ao (and positive) (ao = 0.529 
10 — 10 m is the Bohr radius) it is possible to observe oscillations 
between the usual atomic condensate and the molecular state 
[Tl with a frequency of oscillations of hundreds of KHz Q|. By 
contrast, in the 'Bose-Nova' experiment [TJ typical fields were 
around 167G, the scattering length was only tens of Bohr radii 
(and negative) and the frequency of atom - molecule oscillations 
may be estimated as well over ten MHz (3- While coherent 



The model is based on the Hamiltonian operator for 
N interacting atoms with mass M in a trap potential 
V (r) = (ujIz 2 +ujpp 2 )/2, with radial p and longitudinal 
z coordinates measured in units 2 . of a^o. where a^ 
is a characteristic length of the trap, with associated 
(dimensionless) frequencies to z = u> ax i a i/u> ~ 1/2 and 
uip = uJradiai/^ ~ The interaction is assumed to be 
short ranged. We introduce a dimensionless field oper- 
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ator \& (r) 

constant u = (7Lwa| ) * U = 4?r (a/a^o) 

obeys the equation of motion ^> = 
and satisfies the equal time commutation relations 
[* (*, r) , <f t (t, r')] = (r - r') . We decompose the 
Heisenberg operator \& = $(r,t) + ijj(r,t) into a c- 
number condensate amplitude <f> and a q-number non- 
condensate amplitude ip, consisting of the fluctua- 
tions or excitations. We obtain the equation of mo- 
tion for the fluctuation field by subtracting from the 
full Heisenberg equation the Gross - Pitaievsky equa- 
tion (GPE) governing its own expectation value under 
the self-consistent mean field approximation, ip^ip ~ 



resonance between the atoms and the molecules is expected to 
exist for all of these experiments, and has been shown to play 
an important role in the outcomes of some 0, we deem it un- 
likely that it plays a dominant role in this experiment other than 
renormalizing the scattering length (For details, see |S|). Indeed 
no oscillations are reported in the original experimental paper. 
Instead, as this letter shows, the primary mechanism for the 
Bose Nova phenomena is the parametric amplification of quan- 
tum fluctuations by the condensate dynamics, resulting in bursts 
and jets as particle production from (the squeezing of) the vac- 
uum. Recent numerical simulations |7| and rigorous theoretical 
investigations |S| indicating the inadequacy of mean field theory 
seem to corroborate this view. 

2 We use a sign convention such that the effective coupling 
constant is positive for an attractive interaction, and a sys- 
tem of units where the length a; lo , time ti lo and energy scale 
Eho = ^ w = ^ uj2a ho are defined with reference to the average 
frequency w. We work with units such that these three scales 
take the value 1. 
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(ip^tp) — n, ip 2 ~ (V' 2 ) = m and i/) t ?/; 2 ~ 0. We 
next parametrize the wave functions as $ = <E>oe~ ie , 
ip = -0oe _ie , where $0 and are real. During the 
early stages of evolution, we may regard the conden- 
sate density as time independent, and the condensate 
phase as homogeneous, $0 = $0 (r) , 6 = 6 (t). We 
may then write the equation for the fluctuation field 
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where £0 = 5 (w 2 + 2w p ). To solve this equation 3 we 
decompose it into a self-adjoint and an anti-adjoint 
part ipo = £ + i?7, each part satisfying an equation 
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[H-Eo- 2u*§] £ = 0. 
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Since the trap Hamiltonian is time - independent, we 
have 



dt 2 



+ [H-E Q }H eff ^ = Q. 
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Here H. 



eff 



H — En — 2it$ 2 l - To have an unstable con- 



densate it is necessary that at least one of the eigen- 
values of H e ff is negative; the boundary of stability 
occurs when the lowest eigenvalue is exactly zero. One 
further consideration is that we are interested in the 
part of the fluctuation field which remains orthogonal 
to the condensate, since fluctuations along the conden- 
sate mode may be interpreted as condensate fluctua- 
tions rather than particle loss The ground state 
of H e ff is certainly not orthogonal to the condensate, 
since neither have nodes. 

If we adopt the values u> z = 1/2, u p = y/2, relevant 
to the JILA experiment, then instability occurs when 
k = Nna cr it/ aho = 0.51. This result compares remark- 
ably well with the experimental value k = 0.55 
as well as with the theoretical estimate presented in 
Ref. ^2] • This agreement may be seen as natural, as 
the equations we postulate for the fluctuations may be 
obtained from the linearization of the GPE, discarding 
both n and rh. In both calculations, the geometry of 
the trap plays a fundamental role. 

3 The squeezing of quantum unstable modes and its back re- 
actions on the condensate has been considered before, e.g., as a 
damping mechanism for coherent condensate oscillations Q], but 
the application to the description of condensate collapse has up 
to now been mostly qualitative 1101 . 



Scaling of t co u apse and Critical Dynamics As we 

have already noted, even for condensate densities above 
the stability limit, no particles are seen to be lost 
from the condensate during a waiting time t co u apse . 
Experimentally, t co u apse is seen to get very large when 
the threshold of stability is approached from above, 
in a way which closely resembles the critical slowing 
down near the transition point characteristic of crit- 
ical dynamics. In our problem, the quantity which 
plays the role of relaxation time is the characteristic 
time e _1 of exponential growth for the first unstable 
mode. This quantity diverges at the stability thresh- 
old, which in our analogy corresponds to the critical 
point. By dimensional analysis, we are led to the es- 
timate t co u apse ~ e _1 . Close to the critical point, we 
find 
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The power law Eq. describes with great accuracy 
the way t co u apse scales with the scattering length; the 
best fit to the experimental data is obtained for t cr i t ~ 
5ms. 
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Figure 1: 

In Fig. 1 we plot the scaling law (JjjJ (full line) derived 
here and compare it with the experimental data for 
Nq = 6000 as reported in Refs. (small black points), 
the ijvi ~ (uNo) prediction (suitably scaled) as 
given in |101 IT5| (dashed line) and the results of nu- 
merical simulations reported in j!4| (large grey dots). 
While all three theoretical predictions may be consid- 
ered satisfactory, the t^L ~ {uNo)~~ l behavior fails to 
describe the divergence of t co u apse as the critical point 
is approached, and the results of numerical simulations 
reported in |14| based on an improved Gross-Pitaevskii 
equation tend to be systematically above the experi- 
mental results, which may be a further indication of 



2 



the quantum origin of this phenomenon. 

Bursts and Jets as Amplified Quantum Fluc- 
tuations We now consider the evolution of quantum 
fluctuations, treated as a test field riding on the col- 
lapsing condensate whose dynamics is extracted from 
experiment. The initial state is defined by the condi- 
tion that u = for t < 0; we shall take it to be the 
particle vacuum |0), defined by ipo (x,0) |0) = every- 
where. 

One can introduce a mode decomposition of the £ op- 
erator based on the eigenfunctions of [H — Eo]H e tt. 
For short wavelengths A, since H ~ A -2 >> 2u$q, 
we expect these eigenfunctions will approach the trap 
eigenmodes. The fact that particles in bursts are seen 
to oscillate with the trap frequencies pQ also suggests 
that their dynamics is determined by the trap Hamil- 
tonian. Based on these observations we can assume a 
homogeneous condensate 2it<i> 2 l ~ K^ 1 aut z No (t), where 
Ao (t) is the instantaneous total number of particles 
in the condensate. In practice, k _1 is a measure of 
the overlap between the condensate and the excitation 
modes. Therefore, the approximation may be improved 
by adjusting k according to the range of modes where 
it will be applied. 

Let No be the initial number of particles in the 
condensate, and a cr = k/Nq the corresponding crit- 
ical scattering length. Trap eigenfunctions ipn (r) 
are labelled by a string of quantum numbers n = 
(n z ,n x ,n y ) . The eigenvalues of the trap Hamiltonian 
are (with the zero energy already subtracted) Ef L — 
u z n z + Up (n x + n y ). There are two kinds of modes, 

stable (oscillatory, or thawed) modes if Eft > (^f - ) 
and unstable (growing, or frozen) modes if not. In the 
former case we find that, although we assume vacuum 
initial conditions, these modes do not remain empty. 
Up to tcoiiapse-, when the number of particles in the 
condensate is constant, the density 
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(where = y Eft Eft — (^f - ) Uz ) nas a cons tant 
term and an oscillatory term. This oscillatory term is 
responsible for the appearance of 'bursts' of particles 
oscillating within the trap observed in the Bose-Nova 
experiment 1 j. In the WKB limit it describes a swarm 



of particles moving along classical trajectories in the 
trap potential. 

In the opposite case Eft < (jr~j w zi the formulae 
for the density is obtained by the replacement of 
in © by ian, thus loZ 1 sinu^i — > aZ 1 sinh ovjt. Physi- 
cally their difference is immense. In the first place, the 
density is growing exponentially, but unlike the previ- 
ous case, there is no oscillatory component, and these 
particles do not oscillate in the trap, in the sense de- 
scribed above. These modes come alive at r evo i ve (as 
the scattering length is set to zero), whence they be- 
come ordinary trap modes which oscillate in the trap 
in the same way as the the burst modes . To the ob- 
server, they appear as a new ejection of particles from 
the core of the condensate, which makes up the so- 
called 'jets'. The sudden activation of a frozen mode 
(we are borrowing the language and concept of cosmo- 
logical structure formation) by turning off the particle 
- particle interaction may be described as a "thaw" . 

Observe that in this picture several conspicuous fea- 
tures of jets become obvious. Jets may only appear if 
the turn - off time r evo i ve is earlier than the formation 
time of the remnant. Once the condensate becomes 
stable again, there are no more frozen modes to thaw. 
On the other hand, jets will appear (as observed) for 
Tevoive < t collapse, when the condensate has not yet 
shed any particles. Also jets must be less energetic 
than bursts, since they are composed of lower modes. 

Beyond t co u ap se the number of particles in the con- 
densate, and therefore the instantaneous frequency of 
the excited modes, becomes time dependent. If we con- 
fine ourselves to the early stages of collapse we may 
assume nevertheless that the condensate remains ho- 
mogeneous. Shifting the origin of time to t co u apae for 
simplicity, we write N (t) = N exp (— i/r) (see Fig. 
2). 




Figure 2: 

After expanding in trap eigenmodes we find the two 
kinds of behavior described above. If Eft > ( ^* ) , the 
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mode is always oscillatory. If E n < , the mode 

is frozen at t co ii apse , but thaws when exp (— t/r) ~ 
Ef{d/auj z . During the frozen period, the modes are am- 
plified, but they only contribute to bursts after thaw- 
ing. If the evolution is interrupted while still frozen, 
they appear as a jet. We therefore conclude that the 
number of particles Nj et in a jet at time T evo i ve is essen- 
tially the total number of particles in all frozen modes 
at that time. This is plotted in Fig 3, for iV = 16, 000, 
^radial = HO Hz, Lo aX iai = 42.7 Hz, a = 36a , and 
K = 0.46 , and compared to the corresponding results 
as reported in pQ. 
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Figure 3: 

We see that the agreement is excellent at early times 
(up to about 6ms). For later times, this model overes- 
timates the jet population. This is due to the fact that, 
by not considering the shrinking of the condensate, we 
are overestimating the overlap between the condensate 
and the fluctuations, thus delaying the thaw. It nev- 
ertheless reproduces the overall slope of particle num- 



ber with tv. 



It should also be remembered that 



we are computing the expected number of particles, 
but in the highly squeezed state which results from the 
frozen period, the fluctuations in particle number are 
comparable to the mean number itself. 

In this letter, we have presented a new viewpoint to- 
wards understanding the salient features in the physics 
of controlled collapse of a Bose-Einstein condensate de- 
scribed in the experiment of [Q, i-e, in terms of quan- 
tum vacuum fluctuations parametrically amplified by 
the condensate dynamics. Our way of thinking here is 
influenced by insights from the quantum field theory of 
particle creation and structure formation in cosmologi- 
cal spacetimes as well as theories of spinodal instability 
in phase transitions. One can conceivably design ex- 
periments with BEC dynamics to test out certain basic 
mechanisms and specific features of quantum proceses 



in the early universe, thus opening a new venue for 
performing 'laboratory cosmology'. 
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